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Abstract 
(N 

^ We consider the mutual Renyi information /'") {A, B) = S^ + S)^ — S^^j^ of 

disjoint compact spatial regions A and B in the ground state of a d+l-dimensional 
0^ conformal field theory (CFT), in the limit when the separation r between A and 

1^ B is much greater than their sizes Ra,b- We show that in general I^^-'^AjB) ~ 

^T C]^ Cj^ {RaRb/^"^)", where a is the smallest sum of the scaling dimensions of 

(y-^ operators whose product has the quantum numbers of the vacuum, and the constants 

^T^ ^AB depend only on the shape of the regions and universal data of the CFT. 

.1^ For a free massless scalar field, where a = d — 1, we show that C)^ Ry[^ is 

^ proportional to the capacitance of a thin conducting slab in the shape of A in d+1- 

2 dimensional electrostatics, and give explicit formulae for this when A is the interior 

of a sphere S'^~^ or an ellipsoid. For spherical regions in d = 2 and 3 we obtain 
explicit results for C'"^ for all n and hence for the leading term in the mutual 
information by taking n — )■ 1. We also compute a universal logarithmic correction 
to the area law for the Renyi entropies of a single spherical region for a scalar field 
theory with a small mass. 



1 Introduction 

Since the pioneering paper of Srednicki [Ij there has been increasing interest in under- 
standing and quantifying entanglement in quantum field theories. In that paper it was 
shown that, in a free scalar field theory, the von Neumann entropy Sa = — Trp^logpA of 
the reduced density matrix p^ describing the degrees of freedom inside a spherical region 
A, which measures the entanglement of the degrees of freedom in A with those in its 
complement, is proportional to the area of its boundary. Subsequently this 'area law' was 
show to be generic in space dimensions (i > 2 [2], and this prompted comparisons with 
black hole physics. 

However in 1993 Holzhey et al. [3] showed that in a conformal field theory (CFT) 
in d = 1 the entanglement entropy of an interval of length R^ goes like log Ra with a 
universal coefficient proportional to the central charge c. (Similar logarithms are now 
understood to occur whenever d + 1 is even, but for higher d these are non-leading with 
respect to the area term j[4j.) Subsequently this logarithmic behaviour was observed 
in numerical studies of critical quantum spin chains whose long-distance behaviour is 
believed to be described by a CFT [5]. A more complete analysis of entanglement in 
1+1-dimensional CFTs was given in Refs. [B]. 

More recently [7| these methods, which involve using the so-called replica trick of 
computing Sa as limit as ra — )■ 1 of the Renyi entropies 

Sa = lim S^ where 5*^ = {I — n)^^ logTrp^ 

n— >-l 

have been extended to the computation of the entanglement entropy Saub between two 
disjoint intervals A, B and the rest of the system in a l-|-l-dimensional CFT. It was 
shown that this encodes all the data of the CFT, not only the central charge. Moreover 
the mutual information, given by the limit as n — )■ 1 of the difference of Renyi entropies 

I^-\AB)^sf + sf-S^-:>^ 

has an expansion of the form 

&\A^B) = Y.C'A\{k,})C's\{k,}) i^) " ^' ■ (1) 

{k,} \ ^ 

Here Ra-, Rb measure the lengths of the two intervals, and r is the separation between their 
centres. The sum is over a set of scaling operators of the CFT, labelled by fcj, with scaling 
dimension x^ , for each replica j. The coefficients are universal and encode information 
about the correlation functions of these operators in the plane. For RaRb/'^^ "^ 1 the 
leading term in ([l]) comes from the case when only two of the Xk are non- vanishing and 
correspond to the lowest dimension operator in the CFT. It was possible [7j to continue this 
result analytically in n to compute the leading term in the mutual information I {A, B). 

The one-dimensional case has also been studied numerically in a number of papers [8] . 

The arguments of Ref. [7] were based on a kind of operator product expansion first 
used by Headrick [9]. In this paper we argue that ([!]) also holds for the mutual Renyi 



entropies in higher dimensional CFTs. However in this case the coefficients are much 
more difficult to compute, and we succeeded in obtaining explicit results for the leading 
term only in the simple case of a free massless scalar field theory, and when A and B are 
spheres. In principle our methods work in any (integer) number of dimensions, but the 
results are most simply expressed in d = 2 and 3. 
In general we find 

where, in 3+1 dimensions, 

^3 
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For the mutual information {n = 1) this gives 93=1^ = 0.266, to be compared with a 
numerical result 0.26 due to Shiba ||12J . 
In 2+1 dimensions 



(n) "■ 

92 ^ 
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which leads to 92 = |, to be compared with the numerical result 0.37 

We are also able to compute the form of the corrections to the leading term, which 

should be 0{RaRb / r'^Y'^'^'^^ and also 0{RARB/r^Y^^ ■ In principle the coefficients are 

calculable. Note that these are more important in d = 2 than in d = 3 which may account 

for the above discrepancy. 

Our methods use the conformal invariance of the massless free field theory, which in 

fact implies that for Ra 7^ Rb the actual expansion parameter is 

RaRb 



r^-\{RA-RBY' 

In this form, our results also apply to the mutual entanglement between the interior of 
a sphere of radius Ra and the exterior of a concentric sphere of radius Rb-, in the limit 

when Ra ^ Rb- 

(2) 
For n = 2 and general shapes we show that the coefficients C\^ are given by the 

capacitance of d-dimensional bodies A, B in d + 1-dimensional electrostatics. 

The mutual information for a free scalar field in higher dimensions has been studied 

in only a few papers. Casini and Huerta ^0\ and Shiba [H] showed that it should decay 

as {RaRb / ^"^Y'^ at large separations (and {RaRb/^'^Y for free fermions ^Q\)- This was 

based on the expression for the density matrix in terms of the correlation functions [13j 

which holds for any system with a gaussian wave functional. However the coefficients must 

still be determined numerically for finite separations and then extrapolation. Recently 

this was carried out by Shiba [12J for the case of equal spheres in d = 2 and 3 space 



dimensions, and also rings and shells [12] as well as some other shapes [^. This shows 
that the mutual information is not extensive, that is, is not given by a double integral of 
a kernel over the boundaries or the bulk of A and B. 

Free and interacting fermions at finite chemical potential (Fermi liquids) have been 
studied in [15], and the entanglement of the radiation field with a dielectric medium in 

m- 

The layout of this paper is as follows. In Sec. 2 we first recall the expressions for the 
various Renyi entropies 5*)^ as path integrals over n copies of M"'"*"^ sewn together in a 
particular way along the boundary of the selected region R to form the conifolds Cj^ ■ We 
then consider the general case of the small Ra,b expansion, and show how the coefficients 
C^^ in this expansion are related to one- and two-point functions on Cj^^. 

The rest of the paper is devoted to the special case of a massless free scalar field theory. 
In Sec. 3.1 we show that for n = 2 the coefficients C\^ are related to the electrostatic 
capacitance of each region. In an Appendix we derive explicit formulas for the case where 
A and B are spherical, or more generally ellipsoidal, by generalising a famous method due 
to W. Thomson to general dimension d. In Sec. 3.2 we obtain results for the spherical case 
for general n and d, using conformal invariance. We show that for d odd the coefficients 
are polynomials in n which can then be continued to ra — t- 1 to find the von Neumann 
entropy. For d even the result is not a polynomial but explicit results may nevertheless 
be obtained for n integer as well as the limit n — >■ 1. In Sec. 3.3 we consider the case 
of a single spherical region in a free massive scalar theory, and confirm the existence of 
universal logarithmic terms in the Renyi entropies first predicted in [6]. 

2 General form of the expansion for Ra,b ^ ^• 

Consider a d-l-l-dimensional quantum field theory in a domain D, in its ground state |0). 
In this paper we consider D to be M.'^, although some of the results may be adapted to 
semi-infinite and finite domains and also to finite temperature. 

If X is some subdomain, we suppose that, in the presence of a suitable UV cut-off, 
the Hilbert space can be decomposed as T-Lx ® T^x- The Renyi entropies are Sx = 
{\ — n)'"^ \ogPx\ where 

Pi") = Tr^^p^ with px = Tr^^|0)(0|. 

As explained in Ref. [6], these can be expressed in terms of a path integral on a 
particular conifold as follows. The ground state wave-functional ({^(a;)}|0) is given by 
a path integral in imaginary time on the half-space ]HI_ = M'^ x (— oo,0) from r = — oo 
to r = 0, with the fields constrained to take the values {^(x)} on r = 0. Similarly 
(0|{'?/'(x)}) is given by the path integral on IHI+ from r = to r = oo. Each of these 
should be normalised by Z~^/^, where Z is the partition function on W^^^ . [Strictly 
speaking we should restrict the integration to |r| < T before taking the ratio, then let 
T — )■ oo, and similarly with the thermodynamic limit in space.] 



To get PJ^ we take n copies of H^ labelled by j = 0, ... ,n — 1, and sew EI_ to 
Erl (mod n) along r = for x E X, while we sew EI_ to Wf for x ^ X. This 

gives a conifold C^ with a (i— 1-diniensional submanifold of conical singularities along 
the boundary dX fl {r = 0}. Then 



Px 



(n) _ Z{Cj^ ] 



In general, for d > 1, Z{C^ ) has a leading term going like exp(a„Vol((9X)A''^^) where 
A is an ultraviolet cut-off, and the dimensionless coefficients a„ are non-universal. This 
gives a term in the Renyi entropies proportional to a„Vol(9X)A'^~^ - the famous 'area 
law' term, which is non-universal and therefore theoretically less interesting. However, in 
the case where X consists of two disjoint compact regions (A, B) the quantities 



/(-)(A, B) ^ 4") + sp - s^i, ={n- ir log L5nf;!5l l (2) 



(n) 
Aul 

z{cf)z{cf' 



are free of these non-universal contributions, and should depend only on the geometry of 
A and B and the universal data of the renormalised QFT. 

In general, however computing Z{C^^]g) is difficult, and, even for d = 1, we have 
explicit results only for free field theories [71 110] . However in the limit when the linear 
sizes Ra,b of A and B are much smaller than their separation r, it was shown in Ref. [7] 
that an expansion in increasing (in general, fractional) powers of RaRb/^'^ is possible 
with coefficients which are calculable in principle. We now generalise this argument to 
dimensions d > 1. 

The basic idea is that, from the point of view of an observer far from A oi B, the 
sewing together of the copies along the boundaries of A and B should be expressible as 
a weighted sum of products of local operators $^"' at some conventionally chosen points 

{fAi fs) inside A and 5, where $^ is in the algebra of local operators of the QFT defined 
on jth copy of W^^^. That is, the sewing operation in each region can be thought of as a 
semi-local operator which couples together the n QFTs, but can itself be expressed as a 
sum of a product of local operators in a direct product of the QFTs. We write 

^y-'MjE) _ /v(")v(")\ cw 

-^ - \^A ^B /(M'^+i)" , [^) 

where 

and similarly for B. Here kj label a complete set of operators on the jth copy. The 
prefactor on the rhs is inserted since we expect the leading behaviour as r — > oo to come 
from the term when all the $fc. are the identity operator, and in this limit ^aub ~ ^a^b- 



The main point now is that the coefficients C*r^' i should be universal and therefore 
independent of the other regions and other local operator insertions as long as they are far 
away. Thus we can compute them in the simpler situation when X = A but we consider 
the correlation functions of an arbitrary set of local operators at points r^^^ on C)^ outside 
A: 



1) 




{kj} j—0 , I /|nd+l-,„ 



Note that the rhs decomposes into a sum of products of 2-point functions on each copy 
of M'^"'"-'^, that is we may take j' = j. 

This is valid for any QFT, but in the special case of a CFT we can choose the complete 
set of local operators so that their 2-point functions are orthonormal for all separations: 

where Xk is the scaling dimension of ^k and we have assumed scalar operators for sim- 
plicity. 

Thus we find, taking the limit r*--'-' — )■ ooj (that is, infinity on the the jth copy of M'^^^) 



C4^= hm Ir^^^l^^'^^^iH^kM'^)), 



(n) 
A 



Note that CA.^ oc R^^ ■' by dimensional analysis. 

(n) 



In 1+1 dimensions, when A is an interval of finite length, C)^ may be uniformized 
to C by a conformal mapping and therefore the CA.i are easily computable, at least for 
the ffist few leading terms in the expansion. In higher dimensions this is more difficult. 
Sec. 3 of this paper will be devoted to the simplest case of a free scalar field theory. 



so that 



However, supposing that we have computed the CA.y, then substituting into (34) and 
again using orthonormality of the operators 

{kj} 

4f^ = y Cf,.Cf,. r-' ^^ ^•'=. , (5) 

pin) p{n) Z_^ iftjl i«jl ' ^ '' 

^A ^B {kj} 

where r = [r^ — r^l. Note that in the ratio all terms which contain the UV divergent area 
law pieces cancel. 

If we now arrange the operators $fc in order of increasing dimension Xk-, this gives 
an expansion in increasing powers of [RaRb/i^'^)- The leading term is unity, and comes 
from taking all the $fc to be the identity operator 1 with x = 0. In rf = 1, the leading 



terms with a single Xk- > in general vanish because the one-point functions of primary 

operators in Cj^ are proportional to those in C where they vanish. However this is 
not necessarily the case for d > 1, unless the one-point function vanishes for symmetry 
reasons. The next contribution comes from taking two of the Xk > 0. Note that, unlike 
the case oi d = 1, these do not have to be equal because orthogonality may not hold on 
C^ ■ However, the leading correction terms will come from the smallest two non-zero Xk, 
and, barring degeneracies, these will correspond to the same operator. 

As an example consider the 2+1-dimensional Ising field theory. The leading operators 
are the magnetisation, with Xo- ~ 0.52, and the energy operator with x^ ~ 1.41. The one- 
point function of the magnetisation on C^ vanishes by the Z2 symmetry of the model, 
but there is no reason for the one-point function of the energy operator to vanish. Since, 
however, 2xo- < x^, the leading term in the mutual information will be proportional to 
{RaRb/'t'^Y^'' , with a correction of order {RARs/r'^Y''. Although the above inequality 
holds for many interacting CFTs, counterexamples exist in supersymmetric theories [17] , 
in which case the leading term in the will correspond to the one-point function on Cab 
of an operator with the quantum numbers of the vacuum. 

A second example is free scalar field theory in all d > 1, to be considered in detail in 
the following section. (For d = 1 the field itself is not a local operator and the leading 
corrections then come from exponentials and derivatives of the field [7].) The field has 
dimension x = {d — l)/2 but has vanishing one-point function, once again because of Z2 
symmetry under — )■ —0. However : 0^ : has dimension d — 1 and a non-zero one-point 
function (see below). The leading term in the mutual information is then a combination 
of these two contributions, and goes like [RaRb / f'^Y~^ ■ For a massless Dirac field, which 
has dimension d/2, the power d—lm this expression is replaced by d. These results agree 
with those of Refs. |10l [I2]. 

The corrections to this leading behaviour come from larger values of the exponent 
'^jXkj in (tSJ), and from higher terms on the expansion of the logarithm in (|2|. For a 
free scalar neld, the first type of correction comes from when either four of the $fc are 
taken to be 0, or when these are taken in pairs as :0^:. These all give a contribution 
0{{RaRb / r'^Y^'^~^^) ■ Note that these are more important for smaller values of d. There is 
also a correction coming from taking one of the $fc to be the stress tensor, which always 
has dimension {d + 1), and taking 2 of them to be the current 9^0. These both lead to 
universal corrections O^^RaRb/^'^Y"^^)- For d > 3 they dominate the other corrections. 

3 Free scalar field theory 

In this section we consider the case of a free massless scalar field. The action is propor- 
tional to J{d(f)Yd'^^^x, and we normalise the field so that its 2-point function in M'^"'"^ 

i£] 

(0(a;i)0(2;2)) = Go{xi - X2) = \xi - 0:21"^'^"^^ • 



^We use r to denote points in M'' and x = (r, r) points in 



As discussed in the previous section, we also need : 0^ :, which may be defined by point- 
sphtting as 

: (P^x) := hm U{x + i(5)0(x - ^5) - Go(5)) . 

Its 2-point function in M'^^^, by Wick's theorem, is 

{■.<f)\xi)::<f)\x2):) = 2Go{xi - X2)\ 

so, to conform to our normahsation convention, we should consider $(a;) = 2^^/^ ■.(j)'^{x) :. 
We would like to compute the analogous correlation functions in C^" . As in Ref. [7], 
rather than thinking of a single free field on this conifold, we think of n copies 0j on M'^"'"^, 
coupled by the boundary conditions across r = 0: 

0,(r,O-) = 0,+i(r,O+) (re A); 

= 0,(r,O+) (r^A). 

Thus the coefficients CA.y we need to lowest order are 



C,y = C'o^...,i,...,i,...,o = lim {xiX2Y~'{(j)j{x,)<j)j,{x2))ci.) 



(6) 



where the non-zero entries occur at j ^ j', and 



^n — ^o,...,2,...,o 



2-^/Mima;2(^-i)(:0?(a:):U). 

■r- — i.rv--i -^ ^ A 



(7) 



In the language of c/+l-dimensional electrostatics on Cj" , {(f)j{xi)(f)ji{x2)) is the potential 
at Xi on copy j due to a unit charge at X2 on copy j', while {: (j?-,{x) :) is the excess 
self-energy of a unit charge at x on copy j' . Note that 

^(0j(a;i)0j/(x2))^w + (:0j,(a;):)^(n) =0, 



which follows from conservation of electric fiux. 
The leading correction in (|5| is then 



r 



-2(d- 



''(ki:cfrCf,' + Y.cfff>) 



"n-l 



i^-2(d-l)^ 



.7=1 ^ 



(0j(xi)0o(a;i))r(") 



+ hm (xf)'^-^(:0^(xi):),(„, lim (xf)'^-^(:0^(xi) :)^,„) 



lim (xfY ^(0j(xi)0o(a;i))^(n) 



(8) 



where we have used the cyclic symmetry to extract an overall factor of n. 



3.1 The case n = 2 

For n = 2 it is useful to define the linear combinations 0-|- = 2~^/^(0o =t 0i), which satisfy 

0„(r,O-) = -0„(r,O+) {re A); 
= +0_(r,O+) (r^A), 

while 0+ is continuous across r = 0. Note that 

(0±(x)0±(xi)) = (0o(a;)0o(a;)) ± (0i(a;)(/)o(a;)) , 

so that the correlation functions on the left hand side can be interpreted as the potential at 
X due to a unit charge at Xi. For the upper + sign, the potential is continuous everywhere 
else, and so is equal to Go{x — Xi). 

For the lower — sign, however, it is constrained to change sign across A (1 {t = 0). 
We notice that in (617) we may take xi — )■ oo in any direction. For convenience choose 



it to lie in the hyperplane r = 0. Then the potential due to a unit charge at xi must be 
symmetric under reflection r — t- —t. Therefore the potential on A fl (r = 0) vanishes. 
Thus, as far as 0_ is concerned. An {t = 0} acts like a conductor, at zero electrostatic 
potential. 

Thus we have the electrostatics problem of finding the potential at x due to a unit 
charge at xi, in the presence of a conductor held at zero potential at A fl {r = 0}. In 
general this is complicated, but since we are only interested in the far field in the limit 
when \xi — r^l ^ Rai we can make a simple approximation, valid in this limit. Define 
(j){x) = {(f)_{x)(l)_{xi)) — Go{x—Xi). Then 0(x) is regular at Xi and takes an approximately 
constant value — |xi|~^'^~^'* on the conductor. This will induce a total charge — Cyi|xi|'^~^ 
on the conductor, where C^ is its electrostatic capacitance. Therefore, as x,xi — ?■ cxo, 
(0_(x)0_(xi)) - Goix - xi) ~ -Ca|x|-('^-i)|xi|-('^-i). Thus 

{4>i{xi)4>o{xi)) ~ ^Ca\xi 
{-.(plixi):) = lim 



¥Xi 



-2(d-l) . 

1 








(9) 


o(x)0o(2;i)) - 


- Go{x - 


-^l)) = 




(10) 



givmg 

^'='(AB)~^. (11) 

This is valid for any compact regions A and B. On dimensional grounds Ca,b oc -RJ"^, 
but the coefficient depends on the shape of the regions. Very few cases are known exactly. 

In the Appendix we show, generalising a result of W. Thomson (Lord Kelvin), that 
when dA is a hypersphere 5''^"^of radius Ra , so that A fl {r = 0} is a disc, 

r(rf/2)r(i/2) ,_, 

^^-7rr((d+l)/2)^^ • ^^^^ 

For d+ 1 = 3 this gives Thomson's result -^Ra, while for c? + 1 = 4 it gives ^R\- 



3.2 Free field theory when A and B are spherical, general n. 

The above symmetry argument does not seem to generalise to larger values of n, but 
further analytic progress can be made in the case when A and B are the interior of 
spheres S'^~^. In that case we can exploit the conformal invariance of the free field theory 
to compute the coefficients C •', . 

Before doing this, we note that conformal invariance implies in this case that /*■"•' (A, B) 
is a universal function of the quantity 

RaRb 



r^-\{RA-RBY' 

Given any two spheres Sd-i A and B, we may expand them into spheres 5*^ of the 
same radii in rf + 1 dimensions about their common equatorial plane r = 0. Confor- 
mal transformations in d + 1 dimensions will transform them into other spheres (count- 
ing hyperplanes as spheres through the point at infinity.) Since the system has axial 
symmetry about the line joining their centres, we may restrict to conformal transforma- 
tions which preserve this line. Under such transformations, the cross-ratio of the points 
(— |r — \Ra, ~|^ + \Ra^ \'r — \Rb) \'r + \Rb) where the spheres intersect this line is 
invariant. This is the quantity above. Since by O) /'^"^(A, B) is given in terms of a ratio 
of partition functions in which all metrical factors cancel, it should be both scale and 
conformally invariant. 

We now use conformal invariance to compute correlation functions on C^ when A is 
a sphere. We may regard A as the intersection of a d-dimensional ball of radius Ra with 
the equatorial plane r = 0. Consider the effect of making an inversion in M"'"'"^ which 
sends a point on the boundary of A to the point at infinity. This maps the boundary of 
the ball S*^ into a hyperplane W^ . The plane r = is preserved by the mapping, and 
so the boundary of Ca is mapped into a hyperplane M"^"^. A fl {r = 0} itself is mapped 
into a (i-dimensional half-space. This is easier to visualise in d = 2, when A fl {r = 0} 
is mapped into a half-plane with an infinite line M as its boundary. In the replicated 
theory this turns into a line of conical singularities. The conifold Cj^ is mapped into 
C 2 = {2-dimensional cone of opening angle 2Txn] x M'^^^. 

The correlation functions transform covariantly under this conformal mapping: 



{(t)j{Xl)(f)j,{X2))Ar^) 



dx\ 



dxi 



(d-l)/2 



dx2 



dXn 



{d-l)/2 



{(f)j{x[)(f)j,{x'2)) ^ 



Ca 



) • 



The mapping brings the points at Xi^2 = cxo to a finite distance 1/{2Ra) from the conical 
singularity. The Jacobian cancels the factors of (xiX2)'^~^ in (6|7). Thus 



Cfj, = {(J),{1/2Ra)M^/'^Ra))cM) = {2RAy-'{<l)j{l)(Pf{l))c,M , 

and similarly for Cj^. 

Thus we need to compute the potential on the jth copy at unit distance from the hy- 
perplane of conical singularities due to a unit charge in the same position on the j'th copy. 

9 



As before, because of the cyclic symmetry we can take j' = 0. Since this problem now 
has axial symmetry the calculation is simplified. One approach is to introduce cylindrical 
polar coordinates (p, 6, z), where 6 G [0, 27rn] and z is a {d — l)-dimensional coordinate in 
the W^^^ subspace. We need the Green's function satisfying 

-V2g(")(p, 9, z) oc 5{p - l)5{e)5''-\z) , 

withG'(")(p,^ + 27rn,f) = ^("^(p,^,^). We then have (0j(l)0o(l))c/(") = ^("^(1, 27rj7n,0). 
An expression for G{p^ 6, z) may be found by Fourier transforming with respect to z 
and solving in terms of Bessel functions, but the resultant integrals and sums are ill- 
conditioned and we have not been able to make the continuation in n. 

We adopt a different approach, for which the complete answer for all n may be obtained 
for all d. Instead of considering n to be initially a positive integer, suppose n = 1/m where 
m is a positive integer. The solution for < ^ < 27r/m is then immediate by the method 
of images: 

m—l 

G^^'"'\p, e,z) = Y^ Go{p, e + 27rA;/m, z) . 
Specialising to p = 1, 6* = 0, 

m— 1 .. 

G'(i/'")(1, ^, 0) = y —^ . (13) 

t^(2-2cos(^ + 27rA;/m))('-')/' 

For d — 1 even this sum is straightforward, but more difficult for the odd case, as we 
illustrate below. 

3.2.1 Thecaserf = 3 

In this case the sum can be evaluated explicitly in a number of ways. For example, we 
can regularise it and write it as 



m— 1 ^ m—l oo oo 

hm V = = lim V V V (f+P' 

p^l- /-^ (l - Q(^i(e+2-nk/m)\ \ (l _ ^g-^(9+27rfc/m)^ p^i_ Z-^ Z-^ Z-^ ' 

fc=0 \^ 1-'^ ))\'- H^ ) "^ fc=0p=0p'=0 



^i(p-p'){e+2-Kk/m) 



where \p\ < 1. The sum over k vanishes unless p — p' = (mod m), when it gives m. For 
p > p' we can write p = p' + Im, and similarly for p' > p. Subtracting off the k = term 
to avoid double counting gives 



oo oo oo / -. 1 \ 

V V nV+«mg*«me , p p , Y- 2p' ^ ^ / ^ , ^ \ 

Z^2-.P +C.C. -h^p I _ p2 yi _ p^ime + I _ pe-^me J ■ 



oo oo 

m 

p'=0 1=0 p'=0 



Taking the limit p — )• 1— then gives the simple result for li = 3 



2 — 2 cos mO 
10 



This valid for m a positive integer. However, since it vanishes exponentially fast as 
m — )■ ±2oo, Carlson's theorem ensures that it has a well-defined analytic continuation, to 
other values of m, in particular to m = 1/n. Note however that this does not make sense 



before performing the sum in (13)! 
First, we note that 

(: 0^(1) :),,(.) = lim [j^j^^^^^^ - J^J^o) = ^^ ■ (14) 

As we show in Sec. 3.3, this relates to the coefficient of a universal term in the Renyi 
entropy of a single sphere in the massive theory. 
The first term in (Isl) involves the sum 



n— 1 -. n— 1 ^ 

$:G<">(i,2.j/,,o)^=ijj: 



i=i j 



n^^ (2-2cos(27rj/n))' 



This may be evaluated for n a positive integer by first regulating it as above, and expand- 
ing in powers of p: 



n— 1 oo oo 

')j/n 



EEE(p+i)(p' + i)^'''''^''^^'""'^'^ 

j=l p=0 p'=0 

The sum over j now gives —1 unless p — p' = (mod n), when it gives n — 1. In this case, 
writing p = p' + In, etc., as before, we get 



oo oo 



2 E E(p + i)(p + 1 + ^^)/^''^'" - E(p + 1)'/^'' - E (p + i)(p' + 1)/^^'' 

p=0 1=0 p=0 p,p'=0 

2 1 \ 1 + p"" 2np" 1 



(l-p2)3 (l_p2)2yi_pn (l_p2)2(l_pn)2 (1 - p)4 ' 

Taking the limit p — > 1 gives, after some algebra, 

J = l 

Once again Carlson's theorem assures us of a unique continuation to non-integer values 
of n. 

Putting all the pieces together we find, for rf = 3 



n 



(n2-l)(n2 + ll) (^2-1) 



I^-\A,B) ^ ^^^r-^ ^- -: ' - -' + ^^.^ (2i..)-(2i.. 



2(n - 1) V 720n4 144^4 

n^-1 ( RaRb V 
15n3(ri - 1) V r^ ) 
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\2 



Note that for n = 2 this agrees with (11), using (12). This is a non-trivial check of our 



methods. Taking the derivative at n = 1 we find for the mutual information the leading 
term 



I{A,B) 



4 [RaR 
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AJ^B 



Clearly a similar calculation can be done for any even d — 1. The result will always be 
a rational function of p and p", which will give, on taking the limit p — )■ 1, a polynomial 



m n. 



3.2.2 The case d = 2 



In this case the sum in (13) is 



m—l 



g(^/'")(i,^,o) = ^ 



t^Q (2-2cos(^ + 27rA;/m)) 



1/2 



m—l _ 



Trfc ^ 
m ' 



(Note that in the physical region < ^ < l-n jm the sine is always positive.) Now use the 
integral representation 



1 



°° x^~^ 



dx 



(0 < p < 1) . 



sm vrp TT jg 1 + X 
Inserting this and performing the sum gives 

1 r x(^/2-)-i(i-x) 



G^^''^\i,e,Q) 



27r Jq (1+x)(1 -xV^) 



dx . 



(15) 



The analyticity properties of this expression in 1/m may be inferred by dividing the 
integration region into (0, 1) and (1, oo) and expanding in powers of x (1/x) in each case. 
This gives 



oo oo 



p=0 q=0 



+ 



il+P+i){l+P+i + ^) (-^+P+^)(-|^+P+^ + l) 



In the physical region for 6 this has an accumulation of poles for negative real m, but is 
otherwise analytic. Moreover it grows like \m\ as m — )■ oo except along the negative real 
axis. We infer from this that it has a unique continuation to the whole m plane apart from 
the negative real axis, in particular to m = 1/n where n > 1, which is found by simply 



setting m = 1/n in (15). We remark in passing that as \n\ — )■ oo the sum is 0{n logra) 



consistent with there being a branch cut along the negative real axis. 
We then see that 



(:0o(l)':).,w 



lim — 

e^o 27r 



oo ^(6l/27r)-l / ^ 



1 + x 



X 



1 



X'' 



1 dx 



1 



r.n-1 



27T Jq (1+X)(1-X") 



dx . 
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By substituting x — )■ x ^, we see that integral is twice its value with an upper limit of 1. 
Forn = 2 we get -(l/27r), consistent with (|9 p0p2| . As n -> 1 



::0o(l)':Lw 



(n 



— 1) 1^ logx 



TT 



{n-\\ 



l-x^ 



OO -. 



TT 



The sum in the first term in rtsl) is 



n-l 



n-l 



5:G<-'(1.2.,.or=-i3 5; 



0=1 



i=i 



OO ^oo 



^0 



VT ;r^(2p + i)2 



:i + a;)(l + y)(l-a;")(l-y") 



in -I). 



dxdy 



1 



OO /-OO 



(l-x)(l-l/)(l-(xyr-l) 



(ixdj/ . 



4vr2 7o 7o {l+x){l + y){l-xy){l-x^){l-y^) 
For n = 2 this is l/47r^, consistent with ( 9|l0[l2 ). As n — ;■ 1 we find 

_!i^ r r '2^M dxdy 

This integral may be done by letting y = u/x and first carrying out the x-integral, yielding 

"-1 r(^o,u)''^^^^ 



47r2 7o (1 - uf 

Putting together the pieces, we then find for the leading term in the mutual information 

for ci = 2 

1 I RaRb 



HAB) 



1 1 ».-2 

2 ■ 6 



r-\2RA)i2RB) 



3.3 Universal correction for a single spherical region 

Although the main purpose of this section is to study the mutual information of two 



spherical regions, it is worth noting that the result in ( 14 ) may be used to find a universal 



logarithmic correction to the area law for a single spherical region in d = 3. 

Suppose A is the interior of a sphere S"^^^ of radius R centred at the origin. For 
the time being we keep d general. Let R = (-R, 0, . . .) in the plane r = 0, and invert 
X = (ri, rj_, r) in d + 1 dimensions, sending the point —R to infinity: 

1 
2R' 



(ri + R)^ + ri + T^ 





ri 


+ R 




in 


+ R) 


2 + rl 


+ r2 


in 


+ R) 


T 


+ r2 
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(:0^(x):) (n) 
The jacobian pre-f actor is 



As before, this maps the conifold C^ into C ]^ = {cone in {r[, r')} x M^ ^, and 

, , d-l 



dx 



{■■(P'^ix'):) ,(n). 
'- A 



and 

(:0,2(a;'):)^,,„)=a„(r;Vr'')-('^-^)/^ 



where, from (14), 



a„ = (:0j(l):)^,w = (l-n^)/(12n^). 



After some algebra we find 



AH 



2 \ id-l)/2 



(:0,(x) :).(„) = an . o , _, poio , ,00,0 • (16) 



cr "^[r2 + r2-i?2] 2 + 4^2^2 



Note that the final answer depends only on r"^ = xf + x^ and r, as it should. At large 
\x\ it decays as {x'^)'^'^^^^ so is integrable at infinity ii d > 3. On the other hand, close 
to the conical singularity (r ~ /?, r ~ 0), writing r = R + u with \u\ <^ R, it behaves as 
(u^ + T'^)-(d-^y'^^ so is integrable if rf < 3. 

Therefore for d = 3 the integral of (: 0^(x) :)^(n) over the conifold is logarithmically 
divergent in the massless theory both in the UV and the IR. The former divergence 
necessitates a short- distance cut-off A~^. The IR divergence may be regulated by assuming 
the theory has a small mass m ^ R^^. However, since both divergences are logarithmic, 
to leading order we may still use the expression (16). This gives from large distances 

^l/m 4^2 



, J —d'x ~ an ■ (4i?2)(2vr2) log(l/m) 



and from short distances 



ani^TiR') [ -^!^^a„(47ri?2)(27r)logA. 



'1/A 

Happily the two coefficients agree, so we get 



(: 0^(f ) ■.)Mn)d^rdT ~ na„(87r^i?^) log(A/m) 



-A 



The extra factor of n comes from summing over j . 
The action for the massive theory is 

S=^ f {{d<pf + m'<i?)d^x. 



57r^ 
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(note the factor of (47r^)~^ inserted to conform with our field normahsation convention so 
that the propagator is 1/x^ when m = 0.) 
Then 



which gives a universal logarithmic correction term to the area law for mR ^ 1: 

a4") = -^i?2^2iog(A/m) . 

Since the dependence on A comes from the region near the conical singularity, we may 
conjecture that this is a special case of 

AS!r' = -^^^™Mog(A/m). 

12n 4tt 

valid for any region A with a smooth boundary dA. This agrees with a calculation in [6J 
for the case when A is a half-space (see also [IQl) based on the 1+1-dimensional result: 

'^+-'-A / r^ A\ /l /I 9. . 2n1/2 "• ^ 



Sa^ ~ ^^Area(M) /log(fc2 + m^ 



(2vr) 



after differentiation with respect to m^. 

We emphasise that this universal logarithmic term is present only in the massive 
theory, and is not the same as the term 0(log(Ai?)) discussed in [1] (see also [TO]). 
Although we have carried out our calculation for d = 3, it may be generalised to any 
even d+ 1. For d+ 1 odd, logarithms of the type we have discussed are absent. 

4 Conclusions 

We have presented analytic calculations of the mutual Renyi information /'■"-'(A, B) of two 
disjoint regions in the ground state of a conformal field theory in general space dimension 
d. They are given as an expansion in powers of the ratio RaRb/t^ with powers which are 
integer combinations of the scaling dimensions of the local operators in the theory, with 
universal coefficients, which factorise between A and B. 

Unlike the case of d = 1 we are able to obtain explicit results only for a free scalar field 
theory, although similar methods should work for other free theories, and even interacting 
theories using the e-expansion. For n = 2 we showed that the coefficient of the leading 
term is proportional to the product of the capacitances of A and B in d + 1-dimensional 
electrostatics. For spheres we computed this by generalising an argument of W. Thomson 
to general d. 

When A and B are spheres we can obtain explicit results for the coefficient of the 
leading term for all n and d, although the computations increase in complexity with d. 
For d+1 even the results are polynomials in n, while for d+1 odd they may be expressed 
as Dirichlet series. We remark that from a technical point of view it was simpler first to 

15 



compute the Renyi entropies for n = 1/m where m is a positive integer, and then perform 
the analytic continuation. This approach may be more generally useful. 

It would be interesting to compare our results with those predicted by the holographic 
approach of Ryu and Takayangi |18j, which were hypothesised to extend to the case of 
more than one region in [19] (however, see [20]). We note that our result for n = 2 that 
the mutual Renyi information is given by the product of the capacitances suggests that 
it cannot be interpreted as an extensive quantity. 

Although we have considered only the massless case for the mutual information, our 
results extend straightforwardly to QFTs with a mass scale m, as long as Ra, Rb ^ m~^, 
by replacing massless propagator r^i'^^'^) by its massive version. We also confirmed the 
existence of a universal correction O(-R^m^logm) for the Renyi entropies of a single 
spherical region in the massive theory with mR <^ 1. 
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of his paper and for further correspondence. I thank Pasquale Calabrese for his comments 
on a draft of the present paper, and Mark Srednicki and Erik Tonni for discussions. I also 
thank Slava Rychkov for drawing my attention to the result in Ref. [T7] . 

A Thomson's argument for general d 

Consider electrostatics in D = d + 1 dimensions. What is the capacitance of a hollow 
ellipsoidal shell with semi-axes Oi, . . . , qdI In the limit a^) — )■ this will give the formula 
for a fiat disc in the shape of a (i- dimensional ellipsoid. We need to know what charge 
distribution produces zero field inside. We generalise an argument due to W. Thomson 
[2T] (see |22| for a modern account.) 



^se+'/'-^^3- 




Figure 1: Newton's demonstration that the field at any point in the interior of a spher- 
ically symmetric charge distribution vanishes. [Reproduced from Newton's Principia, 
www.gutenberg.org/files/28233] . 

Start from Newton's observation |23] that for a sphere, when all the aj are equal, 
the field vanishes when the charge distribution is uniformQJ This may be understood as 

■^Of course Newton was considering gravity, but the argument is valid for any r"'' force law in d + 1 
dimensions. 
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follows: consider a point P inside the sphere, and a bi-cone of small solid angle dQ which 
intersects the sphere in two regions of areas proportional to \IH\'^ and \KL\'^ shown in 
Fig. [TJ The contribution to the electric field at P from these two regions is proportional 

to 

\IH\'^ \KL\'^ 

\PH\'^ ~ \PKY 

directed along HP. But this vanishes by similarity of the two triangles. 

Note that this extends to the case of any spherically symmetric charge distribution. 
In particular, we can consider a uniform charge distribution between two shells of radii a 
and a + 5a. 

Now consider the effect of making a uniform shear transformation Xj — )■ x', = XjXj 
where Yl- Xj = 1. This distorts the two spheres into similar ellipsoids with axes Xja and 
Xj(a + 6a). Gauss' law remains true as long as we rescale the electric field components 
Ej — )■ XjEj, so the electric field inside the ellipsoid still vanishes. The bulk charge density 
between the ellipsoids remains uniform. 

Now take the limit 6a — ?• 0. We need to work out the thickness of the shell at a 
point {x'A on the inner shell, which will be proportional to the surface charge density. 
This came from a point {xUXj} on the sphere, and hence its image on the outer shell is 
x'j{l + 6a/ a). To get the thickness we need to form the inner product of Xj{6a/a) with 
the normal at the point x' . The equation of the ellipsoid is 

j 

so the normal vector satisfies ^ ■ Ujdxj = where dxj is any vector satisfying ^ . Xjdxj/X'j = 
0. The solution is to take rij oc Xj/X'j, and so the thickness at Xj is proportional to 



i 



(E,(4A1)) 



TE^iY^^^vm 



This is the surface charge density required to ensure that the field within an ellipsoidal 
shell vanishes. 

In our case we need to take the limit A^ — )■ oo. Rescale xd = bXo where |6| < 1. Then 
the eUipsoid degenerates into 



and the charge density is 



d 2 



A!, 1 1 



oc — oc 



Xd \b\ I _ d ^ 



~^ 
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For the case of a spherical disc when all the aj = a we get a oc l/y/a"^ — r^ for any d. 
In this case the total charge is 

Q oc r /''' dr = a^-' r'\sm Oy-'dO = a^"^ r(rf/2)r(l/2) 



/o V^^^^ Jo " ' 2r((rf+l)/2)- 

On the other hand the potential is (with the same constant of proportionality) 

V= / , , r-, -^ dr = / de = in/2), 

Jo r'^ Va^ - ^ Jo 



JO I 
SO 



r(^/2)r(i/2) ,. 

^ 7rr((rf+l)/2) 
For D = 3 this reduces to Thomson's result C = (2/7r)a. 
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